arXiv:astro-ph/0211224v2 17 Nov 2002 


Is the afterglow of Gamma Ray Burst 021004 unusual? 

Shlomo Dado^, Arnon Dar^’^ and A. De Riijula^ 

ABSTRACT 

The bumpy light-curve of the bright optical afterglow (AG) of gamma ray 
burst (GRB) 021004, its spectral evolution and its radio emission have been 
claimed to be unusual. In the Gannonball Model of GRBs that is not the case. 
The very early AG’s shape is, as for GRB 990123, a direct tracer of the expected 
circumburst density prohle. The unprecedented precision of the data allows for 
the “resolution” of two cannonballs (GBs) in the AG. These two GBs correspond 
to the two pulses in the GRB and to the two wide humps in the AG light curve. 
The smaller wiggles in the AG are, as for GRBs 000301c and 970508, to be 
expected: they trace moderate deviations from a constant density interstellar 
medium. The observed evolution of the optical spectrum is that predicted in the 
GB model. The X-ray and radio emissions of GRB 021004 are also normal. 

Subject headings: gamma rays: bursts 


Introduction 

Instruments aboard HETE II detected GRB 021004 at 4.50432 UT October 2002 (Shi- 
rasaki et ah 2002). Its bright fading optical afterglow was discovered by Fox et al. (2002) 
nine minutes after the burst and was soon followed worldwide with many telescopes, quite 
continuously during the first few days. The unprecedented density and precision of the data 
are reported in Fig. 1. The redshift of the GRB’s host galaxy, z = 2.33, was first correctly 
determined by Ghornock and Filippenko (2002) and conhrmed by Salamanca et al. (2002), 
Mirabal et al. (2002a), Savaglio et al. (2002), Gastro-Tirado et al. (2002) and Matheson 
et al. (2002). The AG’s brightness allowed precise measurements of its temporal decline, 
polarization, spectrum, and spectral evolution. The optical light-curves deviate from the 
smoothly steepening decline observed in most well-studied optical AGs. There are also “un¬ 
expected” features in the spectrum (e.g., Salamanca et al. 2002; Mirabal et al. 2002b) 
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and its evolution (e.g., Matheson et al. 2002). The spectrum of the radio to submillimeter 
wavelengths was also claimed to be “unusual” (Berger et al. 2002a). 

In this letter we show that -in the CB model- the only uncommon feature of GRB 021004 
is the precision of the measurements. Its most “surprising” properties were anticipated, in 
the sense that there are -again, in the CB model- precedents for all of them in past GRBs. 


1. The Cannonball Model of GRBs 

There is mounting evidence that long duration GRBs are produced in the explosions of 
supernovae akin to SN1998bw (Galama et al. 1998), by the ejection of ordinary baryonic 
matter -essentially ionized Hydrogen- in the form of plasmoids or “cannonballs” (CBs), with 
very highly relativistic Lorentz factors (7 ~ 10^) (Dar and De Riijula 2000, 2001, Dado et 
al. 2002a,b,c,d) but otherwise similar to the ones observed in quasars (Marscher et al. 2002) 
and microquasars (e.g., Mirabel and Rodriguez 1994; Belloni et al. 1997; Rodriguez and 
Mirabel 1999; Mirabel and Rodriguez 1999; Corbel et al. 2002, and references therein). 

In the CB model (Dar and De Riijula 2000, 2001, reviewed in De Riijula 2002a,b) long 
duration GRBs and their AGs are produced in core collapse SNe by jets of highly relativistic 
CBs that pierce through the SN shell and the SN progenitor’s “wind” ejecta. A CB is 
emitted, as observed in /i-quasars, when part of an accretion disk falls abruptly onto the 
newly-born compact central object. Crossing the circumburst shells with a large Lorentz 
factor 7 , the surface of a CB is collisionally heated to keV temperatures and the radiation it 
emits when it reaches the transparent outskirts of the shells —boosted and collimated by the 
CB’s motion— is a single 7 -ray pulse in a GRB. The cadence of pulses reflects the chaotic 
accretion and is not predictable, but the individual-pulse temporal and spectral properties 
are. A long list of general properties (Dar and De Riijula 2001) of GRB pulses is reproduced 
in the CB-model, in which, unlike in the standard “hreball” models, the GRBs’ 7 ’s have a 
thermal-bremsstrahlung —as opposed to synchrotron— origin^. 

A CB exiting the circumburst shells soon becomes transparent to its own enclosed 
radiation. At that point it is still expanding and cooling adiabatically and by bremsstrahlung. 
The hard bremsstrahlung spectrum dominates the very early X-ray AG (for a few tens of 
seconds) with a fluence of predictable magnitude decreasing with time as t“^. All X-ray AGs 
are compatible in magnitude and shape with this prediction (Dado et al. 2002a). 


•^In this (Ghirlanda et al. 2002), as in the relevance of the viewing angle, or the explanation of AG 
features by density inhomogeneities, the “standard” model is phagocytizing the GB-model’s views. 
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After the first tens of seconds, a CB’s emissivity is dominated by synchrotron emission 
from the electrons that penetrate in it from the interstellar medinm (ISM). Integrated over 
freqnency, this synchrotron emissivity is proportional to the energy-deposition rate of the 
ISM electrons in the CB^. These electrons are Fermi accelerated in the CB’s tangled magnetic 
maze to a broken power-law energy distribntion with a “bend” energy eqnal to their incident 
energy in the CBs’ rest frame. In that frame, the electrons’ synchrotron emission (prior to 
attennation corrections) has an approximate spectral energy density (Dado et al. 2002b): 


(p - 2) np me c^ [z//z/b] 
(p - 1) z/b a/1 


where fo is a normalization constant, p ~ 2.2 is the spectral index of the Fermi accelerated 
electrons prior to the inclusion of radiation losses, Up is the ISM baryon density, 7 (t) = 
1 /a/1 — (with P = v/c) is the Lorentz factor of the CBs, and z/b is the “injection bend” 
frequency in the CB’s rest frame®. 


z/b ~ 1.9 X 10® /(t)]® 


Ur 


10-3 


cm 


-3 


1/2 
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( 2 ) 


The radiation emitted by a CB is Doppler-shifted and forward-collimated by its highly 
relativistic motion, and redshifted by the cosmological expansion. A distant observer sees a 
spectral energy flux: 




CB 


V, t'], 


( 3 ) 


where R is the radius of the CB (which in the CB model tends to a calculable constant value 
Rmax of C(103"^) cm, in minutes of observer’s time), A(z/, t) is the total extinction along the 
line of sight to the CB, Dl(z) is the luminosity distance (we adopt Hq = 65 km/(s Mpc), 
Dm = 0.3 and Da = 0.7), z/' = (1 + z) z//5(t), t' = (5(t) t/(l-|-z), and (5(t) is the Doppler factor 
of the light emitted by the CB. In the domain of interest for GRBs, 7 ^ 3> 1 and 1, 

with 9 the angle between the CB’s direction of motion and the line of sight to the observer. 


(5(t) ~ 


2 7(t) 

l + 02^(t)2 ■ 


( 4 ) 


"‘The kinetic energy of a CB is mainly lost to the ISM protons it scatters; only a fraction < me/mp is 
re-emitted by the incorporated electrons, as the AG. 

^This bend frequency is not the conventional “cooling break”. It reflects an injection bend in the electron 
spectrum at the energy Eb = 7 (t) me c^ with which the ISM electrons enter the CB at a particular time in 
its decelerated motion (Dado et al. 2002b). 
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The total AG is the sum over CBs (or large individual GRB pulses) of the flux of Eq. (3). 


For an ISM of constant baryon density np, the deceleration of a GB results in a Lorentz 
factor, 7 (t), that is given by (Dado et ah 2002a): 


7 = 


G = 


B = 


7(7o,0,Xoo;t) = - 


e^ + ce^ + - 

1/3 


LB2 + 206 + BVB2 + 406 
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+ 


3 02 
7o 


+ 


Oct 


(5) 


(1 + Z) Xoo ’ 

where 70 = 7 ( 0 ), and Xoo = N(.g/(7r%) characterizes the GB’s slow-down in terms of 
its baryon number and Rmax, hs radius (it takes a distance Xoo/ 70 , typically of 0(1) 
kpc, for the GB to slow down to half its original Lorentz factor). 


The extinction, A(z/,t) in Eq. (3), can be estimated from the difference between the 
observed spectral index at very early time when the CBs are near the SN and that expected 
in the absence of extinction (Fobs oc for z/ z/b). The GB model predicts —and 

the data conhrm— the gradual evolution of the effective optical spectral index towards the 
constant value ~ — 1.1 observed in all “late” AGs (Dado et al. 2002a). The “late” index is 
independent of the attenuation in the host galaxy, since at t > 1 (observer’s) days after the 
explosion, the GBs are typically moving in the optically-transparent halo of the host galaxy. 


The comparison of the predictions of Eq. (3) with the observations of optical. X-ray 
and radio light-curves and spectra is discussed in Dado et al. (2002a,b). The results —for 
all GRBs of known redshift— are very simple, satisfactory and parameter-thrifty. The GB- 
model results concerning X-ray lines in GRB AGs are also exceptionally predictive, simple 
and encouraging (Dado et al. 2002e). 


2. GRB 021004 in the CB model 

Five properties of AGs in the GB model are particularly relevant to GRB 021004: 

a. For most GRB optical AGs, good fits are obtained by approximating the ISM-density by 
a constant (Dado et al. 2000a,b,c,d). One exception is the optical AG of GRBs “caught” at 
very early times, when the GBs are still moving in a wind-generated density prohle oc l/r^. 
In the GB model the AGs’ flux is proportional to the instantaneous swept-in electron density, 
implying Fobs oc 1/C at very early times, as observed in GRB 990123 and GRB 991216. 

b. The achromatic “bumps” in the AGs of GRB 970508 and GRB 000301c are explained 
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(Dado et al. 2002a,b) by inhomogeneities in the density along the CBs’ trajectory, better 
than by gravitational lensing (Garnavich, et al. 2000). Density changes are expected within 
star formation regions and upon exit from the superbubbles where most of the SNe take 
place (Dado et al. 2000a). 

c. The data on past AGs were course enough or started late enough for the contributions 
of different GBs (which can often be resolved as individual pulses in the GRB phase) to 
coalesce into an AG contribution describable as a single GB or a collection of similar ones. 
But individual GBs have different properties, and may even be emitted at somewhat different 
angles, due to precession of the accretion disk relative to the rotation axis of the compact 
object (Fargion and Salis 1995), as observed in the microquasar SS 433 (Margon 1984). 

d. The spectral shape of the AGs and their time-evolution, in particular their steepening 

at the time-varying frequency z/b(t) of Eq.(2) —towards at late time— is very well 

supported by the data (Dado et al. 2002b). 

e. The excess polarization of AGs above that induced by the ISM in our Galaxy may be 
largely due to the host galaxy’s ISM. In that case, it should be correlated with the extinction 
in the host and decline with time as the GBs move into the halo. 

The 7 -ray light curve of GRB 021004 shows two prominent pulses separated by ~ 30 
seconds (Shirasaki et al. 2002 and http://space.mit.edu/HETE/Bursts/GRB021004/). In 
the GB model, these correspond to two dominant GBs. The good quality of the optical data 
for GRB 021004 —and its double-humped 7 burst— make irresistible the temptation to £t 
its broad-band AG light curves with two CBs, emitted in the same direction 6 relative to 
the observer®, but with otherwise free parameters (normalization, 70 and Xqo). We fix the 
spectral index p in Eq. (1) to the expected p = 2.2, assume a density profile of the form 
n(r) = Uoo [(ro/r)^ -|- 1], and fit simultaneously all the reported NIR, optical and radio data. 

NIR-Optical AG: In Fig. 1 we compare the results of the GB-model’s fit to the data for the 
NIR-optical lightcurves. We have corrected for selective Galactic extinction in the direction 
of GRB 021004 (Schlegel et al. 1998): E(B —V) = 0.06 (attenuation magnitudes Au = 0.33, 
Ag = 0.26, Ay = 0.20, Afj^ = 0.16, Aj = 0.12, Aj = 0.08, Ak = 0.04, and Ajj = 0.01, in 
the various bands), and for absorption in the intergalactic space and in the host galaxy in 
the B and U bands, Ab = 0.25 and Ay = 0.50 mag, estimated from the reported spectra of 
GRB 021004 (Moller et al. 2002; Matheson et al. 2002; Pandey et al. 2002; Bersier et al. 
2002; Holland, et al. 2002; Schaefer et al. 2002). The contribution from the host galaxy to 


®We have also made fits with two different values of the emission angles 0. They do not improve 
signihcantly, and yield similar parameters: two different emission angles are not really necessary. 
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the late AG, normalized to 0.66 ± 0.33 /i-Jansky in the R band (Holland et al. 2002), was 
subtracted. The results are very satisfactory though, as behts a rough approximation to a 
no-doubt very complex system, it is not perfect = 3.5 per d.o.f. for 465 data points, if 
all reported statistical errors smaller than 0.05 mag are hxed to 0.05 mag, to account for 
uncertainties in calibration and in the conversion from magnitudes to Janskies). The broad¬ 
band htted parameters are, 9 = 1.47 mrad, 7o[l] = 1403, 7o[2] = 1259 (implying (5o[l] = 542, 
Jo [2] = 576), Xoo[l] = 25 kpc and Xoo[2] = 620 kpc. To demonstrate the real quality of the 
fit, we have blown up the R-band results in Fig. 2. In the region between ~ 0.5 and ~ 5 
days, the data “wiggle” by as much as 20% around the theoretical curve. It would be easy 
to correct for this by assuming similar deviations of the ISM density relative to the constant 
value adopted at large distances, clearly a moot exercise. 

Optical Polarization: The observed polarization of the optical AG of GRB 021004, which 
is consistent with that induced by the ISM in our Galaxy (Govino et al. 2002a,b; Rol et al. 
2002; Wang et al. 2002), implies that the polarization induced by the host’s ISM is rather 
small and consequently the attenuation of the AG in the host must be small, consistent with 
the GB model £t to the AG. 

X-ray AG: The Ghandra observations by Sako et al. (2002) in the 2-10 keV domain, from 
0.867 to 1.874 days after burst, result in a spectrum dn/dE oc and an AG decline 

dn/dt oc both consistent with the optical observations and with the GB-model’s 

expectations: dn/dE oc E“^-^ and, except at earlier times, an achromatic light-curve. In 
this GRB, as in most others, the temporal and spectral behaviours are similar in the optical 
and X-ray bands, as shown in Fig. 1, but synchrotron emission alone underestimates the 
X-ray fluence (in this case, the optical to X-ray effective index is 0.2 units less steep than 
the separate indices). The GB-model explanation is that in the X-ray domain not only 
synchrotron radiation, but also Gompton up-scattering and line-emission contribute to the 
flux (Dado et al. 2002a,b). There is nothing exceptional in the X-ray data for GRB 021004. 

Radio and Submillimeter AG: In Fig. 3 we compare the results of the GB-model’s 
broad-band fit to the radio and submillimeter data. Unlike the X-ray and optical AGs, the 
radio AGs are sensitive to self-absorption in the GBs, parametrized by a single absorption 
frequency (Dado et al. 2002b). The fit value for the GB that dominates the “late” radio 
AG is z/a[2] = 0.98 GHz, similar to those of other GRBs (Dado et al. 2002b). The fit in 
Fig. 3 is very satisfactory, in particular in view of expected scintillations in the radio data. 
The GB model is seen to correctly predict the early temporal increase (to be followed by a 
later turn-over due to self absorption) and the spectral behaviour (F^, ~ ^-i.i3±o.i5 between 
1.43 GHz and 86 GHz, observed around 5 days after the burst) again akin to those of other 
GRBs: 991208 (Galama et al. 2000), 000926 (Harrison et al. 2001) and 000301c (Berger et 
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al. 2000). GRBs 021004 and 000301c also have very similar redshifts and fluences, dim host 
galaxies, bumpy light-curves and analogous “late” optical AGs (GRB 000301c was optically 
detected after ~2 days by Fynbo et al. 2000). 

The GRB proper. In the GB model (Dar and De Riijula 2000), the rest-frame fluence 
of a GB viewed at a small 9 is amplihed by a huge factor, (5g, due to the Doppler boost 
and relativistic beaming of the radiation (Shaviv and Dar 1995). In Dado et al. (2002a) we 
deduced that the GRB photons of the GRBs of known redshift correspond to a total energy 
release in the GBs’ rest system that is in a surprisingly narrow intervaF, 10 ^"^=*=°-^ erg; and that 
the spread in the “equivalent spherical energies” E^, around their mean, E^, 4 x 10^^ erg, 

is mainly due to the spread of their Doppler factors (deduced from the hts to their AGs) 
around their mean: 5o ~ 1100. For GRB 021004, (5o ~ 550 and the GB-model expectation 
is Eg, Eg(5o/(fo)^ ~ 5 X 10^^ erg, in agreement with the observed ~ 4.8 x 10^^ erg (we 
have corrected the value of Lamb et al. (2002) to a redshift z = 2.335, from their adopted 
z = 1.6). The GRB of GRB 021004 is also entirely normal. 

Spectral Variability: The energy dependence of the injection bend implies that the optical 
spectra, Fobs[^, t]) typically steepen in the first few days, as the bend frequency in the ob¬ 
server’s frame zZb^®(t) = z/b(t) (5(t)/(l-|-z) “crosses” the optical band (Dado et al. 2002b). Let 
Ag^®(t) = c/z/()'^®(t). Given the parameters of our £t to the optical AG, we can predict, by use 
of Eqs.(l), the spectral energy density Fa as a function of the ISM density Up. The two GBs 
of GRB 021004 have similar 70 and (5o but different Xqo: their time evolution and injection 
bends are different. The AG is dominated before and after t ~ 0.5 days by one or the other 
GB, whose properties determine the “early” or “late” spectral evolution. The predictions for 
the spectral ratio at the “late” times (ti = 0.756 d, t 2 = 2.786 d) measured by Matheson et 
al. (2002), are shown in Fig. 4. The upper panel is for no^ = 0.0018cm“^, the density pre¬ 
dicted from Xoo[2] using the GB-model reference values (N^^^ = 6 x 10^°, Rmax = 2 x 10^^ cm; 
Dado et al. 2000 a,b). For this Uoo, Ag'’’®(ti) = 3140 A and Ag’^®(t 2 ) = 8940 A. 

The interpolation between the limiting behaviours at z/<^z/b and z/3>z/b could be more abrupt 
than in Eq. (1). In the lower panel of Fig. 4 we illustrate the result for a sharp transition, 
for which, with Uoo = 0.0048 cm“^, Ag^®(ti) = 1932 A (below the measured range) and 
-^^^( 12 ) = 5500 A, at which point Fa changes from ~ [A/Ab]“^'^ to ~ [A/Ab]“°'®, so that the 
ratio shown in the figure is [A/Ab(t 2 )]“° ® below 5500 A, and unity above. 


^GRBs in the CB model are much better standard candles than in the standard model (Frail et al. 2001). 



3. Conclusions 


We have shown that the CB-model’s interpretation of all the data concerning GRB 
021004 is quite simple and very successful. The only novelty is that, relative to other GRB 
AGs, the optical data start very early and are so copious that a single GB would not fare well 
in explaining the light-curve shapes. Two large pulses are observed in the GRB’s 7 -ray light 
curve, enticing one to fit the rest of the data with two GB contributions. The results are 
excellent and do not require new ad-hoc assumptions about the circumburst ISM density: 
GRB 021004 fits well and naturally along all other GRBs of known redshift. 

Having described a relatively simple AG evolution with 8 parameters, we cannot claim 
that they are all uncorrelated and well determined, although their use results in the impres¬ 
sive prediction for the spectral evolution reported in Fig. 4. We regard our results as an 
existence proof of a very simple GB-model understanding of the bumpy AG of GRB 021004. 

It is, as usual, instructive to compare the GB-model’s results to those of the standard 
“hreball” scenarios. Two of the most striking differences are their degree of simplicity and 
predictive power; the explanation of the alleged X-ray lines in GRB afterglows is an extreme 
case®. The GRB under discussion is also in an extreme class: its standard explanations 
(Lazzati et al. 2002; Nakar et al. 2002; Heyl and Perna 2002) require model-dependent 
inversions of the observed light curves to obtain the input functions (the ISM-density or 
energy-supply prohles and the extinction in the host) needed to “predict” the output func¬ 
tions (the light curves at different frequencies). The various quoted authors extract very 
different input functions, thereby illustrating the limited “posdictive” power of hreballs. It 
would be difficult, on these grounds, to falsify any model with that much freedom. 

For a non-constant ISM density, short-time spectral variations and AG features are 
expected in the GB model due to the local instantaneous dependence of F,^ and z/b on the 
density, Eqs. (1,2). These are not expected in the standard models because signals emitted at 
a hxed t from different portions of the fireball or hrecone are observed at different times and 
with different Doppler boosts, erasing all effects of width At t. Moreover, the standard 
AG fiuence depends on the integrated density (up to time t): an abrupt rise of fiuence such 
as that observed in GRB 970508 requires a different explanation (the AG is “re-energized”, 
as in Piro et al. 1998). A similarly abrupt dip would be unexplainable. 

The GB model is very modest in the adjectives that refer to GRBs. None of them 
is exceptional, not even the very energetic GRB 990123, nor 970508 with its peculiar AG 


®In the CB model the existence and time-dependent energies of the lines are predicted, requiring no 
custom-made metallic envelopes, nor any other ad-hoc assumption (Dado et al. 2002e). 
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shape, nor the extraordinarily close-by 980425. They are all associated with supernovae, 
seen when they are visible, not seen when they are not (Dado et al. 2002a). The explosions 
that generate GRBs are not “the biggest after the Big Bang”. The mechanism that begets 
GRBs is common; it takes place in quasars and microquasars as well. The model works very 
well and is very predictive, thus falsihable. 

Acknowledgment: This research was supported in part by the Helen Asher Space 
Research Fund and by the VPR fund for research at the Technion. Useful remarks by 
Stephen Holland are gratefully acknowledged. 
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Fig. 1.— Comparison between the NIR-optical observations of the AG of GRB 021004 and 
the GB model £t for two GBs with different parameters, each contributing as in Eqs. (1) to 
(5) and corrected for Galactic, intergalactic and host extinction as described in the text. The 
ISM density was assumed to be a constant plus an additional “wind” contribution decreasing 
as 1/r^. The hgure shows (from top to bottom) 10^ times the K band, 10^ times the H-band, 
10^ times the J-band, 10 times the I-band, the R-band, 10“^ of the V-band, 10“^ of the B- 
band and 3 x 10“^ of the U-band. The observational data points are those reported to date, 
in GGN notices (recalibrated with the observations of Henden et ah 2002), and in Bersier et 
ah (2002), Holland et ah (2002) and Randy et ah (2002). An estimated contribution of 0.66 
p-Jansky of the host Galaxy, (Holland et ah 2002) was subtracted from the late-time data 
in the I, R and V bands. The bottom line is the GB-model prediction for the synchrotron 
contribution to the 2-10 keV X-ray AG (Datum is from Sako et ah 2002). 
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DAYS AFTER BURST 


Fig. 2.— Blow-up of the R-band results of Fig. 1. The ISM density was assumed to be a 
constant plus an additional “wind” contribution decreasing as 1/r^. The wind contribution 
is only signihcant at t < 0.01 days, after which the CBs are more than 1 pc away from the 
progenitor. The “wind” is also seen in other AGs observed early enough (Dado et ah 2002a). 
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DAYS AFTER BURST 


Fig. 3.— Comparison between the radio and submillimeter observations of the AG of GRB 
021004 and the GB-model light-curves calculated from the broad band £t to its AG with 
two GBs with different parameters, each contributing as in Dado et ah (2002b). The hgure 
shows (from bottom to top) the data and the predicted lightcurves at 1.43, 4.86, 8.46, 15, 
22.5, 86.23 and 231.7 GHz. Data points are from Frail et al. (2002), Berger et ah (2002a,b); 
Pooley et al. (2002a,b,c) and Bremer et al. (2002). The symbols for the data are the same 
as those for the corresponding theoretical curves. 
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Fig. 4.— Comparison between the ratio of the spectra (F^) of the optical afterglow of GRB 
021004 at ti = 0.756 days and t 2 = 2.786 days after burst, normalized to ~1 at long wave¬ 
lengths, as measured by Matheson et ah (2002). The wiggles on the curve extracted from 
observations are artifacts of the smoothing procedure, the error in the plotted magnitude 
difference is estimated at 0.1 mag. The smooth line in the upper panel is the prediction 
from Eqs. (1,2), with the ISM density extracted from the observed Xoo[2]. The smooth line 
in the lower panel illustrates an abrupt injection bend at Ab(t 2 ) = 5500 A. 






